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Abstract 

Stationary periodic solutions of the two-dimensional Gross-Pitaevskii equation are 
obtained and analyzed for different parameter values in the context of the problem 
of a supersonic flow of a Bose-Einstein condensate past an obstacle. The asymptotic 
connections with the corresponding periodic solutions of the Korteweg-de Vries and 
nonlinear Schrodinger equations are studied and typical spatial wave distributions are 
discussed. 



1 Introduction 

The Gross-Pitaevskii (GP) equation plays a prominent role in the description of nonlinear 
dynamics of Bose-Einstein condensates (BEC) (see, e.g., PQ). It describes, in the so-called 
mean-field approximation, the behavior of the order parameter ip(r) (the "condensate wave 
function") and has the form 

ih^ = -^-A4, + V{v)4, + g\^ 1 (1) 

where V(r) denotes the potential of the external forces acting on the condensate (e.g. the 
confining trap potential), g is an effective coupling constant arising due to inter-atomic 
collisions with the s-wave scattering length a s (positive for repulsive interactions), g = 
4:7ih 2 a s /m, m being the atomic mass. The GP equation (JTJ takes into account the dispersive 
and nonlinear properties of the condensate which can give rise to various nonlinear structures 
in a BEC flow. In particular, vortices, bright (for g < 0) and dark (for g > 0) solitons have 
been extensively studied theoretically in the framework of the GP equation and observed in 
experiments (see, e.g., [T] and references therein). 

An important insight into the structure of the solutions of the GP equation is provided 
by the fact that its potential-free, one-dimensional reduction coincides with the integrable 



nonlinear Schrodinger equation 

#r + 4>XX ~ 2^1 V = • (2) 

Here X = x/£, r = tc s /£ and o = sgng, where £ = fij ^2mn Q g is the healing length and 
c s = h/\/2m^ is the sound velocity in a BEC of density n . 

In recent experiments j2] on free expansion of a non-rotating BEC after its release from 
a trap, the new interesting blast wave patterns have been observed. Using the modulation 
solutions of the one-dimensional NLS equation (JSJ) (defocusing case), obtained earlier in 
|U E] , and numerical simulations in 2D and 3D, these blast waves have been identified 
in [3 |Hl E] with expanding dispersive shock waves, which represent oscillatory counterparts 
of classical evolutionary gas-dynamic shocks. 

Yet another type of nonlinear wave patterns has been observed in another series of ex- 
periments on the flow of a non- rotating BEC past macroscopic obstacles also reported in . 
In these structures have been associated with spatial dispersive shock waves. Spatial 
dispersive shock waves represent dispersive analogs of the the well-known viscous spatial 
shocks (oblique jumps of compression) occurring in supersonic flows of compressible fluids 
past obstacles. In a viscous fluid, the shock can be represented as a narrow region within 
which strong dissipation processes take place and the thermodynamic parameters of the 
flow undergo sharp change. On the contrary, if viscosity is negligibly small compared with 
dispersion effects, the shock discontinuity resolves into an expanding in space oscillatory 
structure which transforms gradually, as the distance from the obstacle increases, into a 
"fan" of stationary solitons. The theory of the generation of such spatial dispersive shock 
waves has been developed in JJJ in quite general terms for supersonic dissipationless flows 
past a slender body when the flow can be asymptotically described by the Korteweg-de Vries 
(KdV) equation and an effective description of dispersive shocks becomes possible via the 
Whitham modulation theory [T2*l IT3| IT3j. 

In a different approximation, for highly supersonic flow past slender body, for a certain 
range of parameters, the GP equation asymptotically reduces to the one-dimensional non- 
linear Schrodinger (NLS) equation (J2J), albeit for a completely different set of independent 
variables [TO] - As a result, the problem of stationary dispersive shock waves in this case can 
be treated by similar methods of Whitham's theory. It is essential that in all cited papers, 
the analytical advances in the description of the dispersive shock waves have become pos- 
sible owing to a complete integrability of the KdV and ID NLS equations resulting in the 
possibility to represent the associated modulation Whitham systems in Riemann diagonal 
form (see [T3| and references therein), which dramatically simplifies further analysis. 

However, in real experiments the obstacles cannot be treated as slender bodies and the 
flow is not highly supersonic, so the spatial dispersive shock waves must be considered in 
the frame of the full, unapproximated GP equation (|T|I. In such a study, one inevitably 
faces the major obstruction to obtaining analytical solutions to an initial or boundary value 
problem because the multi-dimensional GP equation, even in the absence of the external 
forcing term, is a non-integrable system and powerful spectral methods (inverse scattering 
transform, finite-gap integration) are not available for it. However, the Whitham modulation 
approach still remains a possibility provided a minimal structure (availability of a certain 
number of conservation laws and a traveling periodic wave solution) is present |16j . 
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In the framework of the Whitham modulation theory, a dispersive shock wave is con- 
sidered as a modulated nonlinear periodic wave whose parameters change slowly on a scale 
about one wavelength and one period (see JT3JJT5] for instance). For a certain spatio-temporal 
domain, the dispersive shock represents a train of solitons well separated from each other. In 
the context of the BEC flow past an obstacle a 2D spatial dispersive shock wave is asymptot- 
ically (far enough from the obstacle) represented as a "fan" of spatial dark solitons. When 
the obstacle is not very large, only one soliton can be generated. This simplest case has 
been considered recently in [T]\ where exact spatial soliton solution of the GP equation (JTJ) 
was found for condensate with a flow and it was shown by a numerical simulation that such 
"oblique" solitons can be generated by a supersonic flow past an obstacle. To extend this 
theory to "multi-soliton" dispersive shocks represented by a modulated periodic wave, it is 
necessary, first of all, to find stationary periodic solutions of the GP equation (0), and this 
is the main aim of the present paper. We also note that the periodic solutions of the GP 
equation could be important in other, than BEC, areas where the multi-dimensional NLS 
equation is used to model nonlinear wave propagation. 



2 Periodic solution 

The stationary solutions of Eq. (£Q) can be sought in the form 

V>(r) = v 7 ^) exp U J u(r')dr'j exp (-jr*) , (3) 

where n(r) is the density of atoms in the BEC, u(r) denotes its velocity field and \x is the 
chemical potential. It is convenient to introduce the dimensionless variables 

f = r/v / 2£, n = n/n , u = u/c s . (4) 

Substituting Eq. (J3J) into (JTJ) and separating real and imaginary parts we obtain a sys- 
tem of equations for the density n(x, y) and the two components of the velocity field 

u = (u(x,y),v(x,y)), 



(nu) x + (nv) y = 



n l + n l n xx + n, 



i j, y • "xx i ,u yy i n 

uu x + vu v + n x + - 2£ = 



2 



uv x + vv v + n v + 



An J x ' (5) 



nl + nl n xx + n 



!! ' ;/ \ 8n 2 An 

where we have omitted tildes for convenience of notation. 

One can see that the (0 + 2) reduction (jSJ) of the GP equation (JTJ) is drastically different 
from its (1+1) NLS reduction (J2J) (which also can be represented in a hydrodynamic form 
by the change of variables analogous to Q - see Section 4). Putting aside the subtle 
integrability aspects (the equation (J2J) is a completely integrable system while there is no 
indication of integrability for the system (JHJ)) we note that, first of all, the scalar system 
corresponding to the one-dimensional NLS equation (J2J) consists of two equations (for n and 
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u) while the spatial (0+2) case leads to three equations (jSJ). As a result, the structure of the 
system (j5J) is significantly more complicated. This can already be seen by comparing the 
dispersionless limits of (J2J) and Indeed, the dispersionless limit one-dimensional NLS 
equation (J2J) coincides with the classical shallow-water equations (or classical gas-dynamic 
equations with the adiabatic index 7 = 2) which are always hyperbolic. Contrastingly, the 
dispersionless limit of (jSj) coincides with stationary two-dimensional gas-dynamic (7 = 2) 
equations which have a mixed elliptic-hyperbolic structure depending on the absolute value 
of the flow velocity. Thus, one can expect a considerable difference in the behavior of the 
periodic solutions and their modulations in these two different reductions of the GP equation. 
We look for the solution of the system (jSJ) in the form of a "travelling" wave 

n = n(6), u = u(6), v = v{6), (6) 

where 9 = x — ay, a being the "slope" parameter of the stationary wave (the wave crests lie 
on parallel lines with the slope a to y-axis). Under this ansatz, the first equation (jHJ) gives 
at once 

u — av = —, (7) 

n 

where A is the integration constant, and the other two equations reduce to 

n\ - 2nn K + 2n 3 - j^n 2 + = 0, (8) 

s 1 + a 1 1 + a 1 

where B is another integration constant and we have introduced new independent variable 

26 2(x-ay) 
Z = , = . (9) 

VTTa7 VTTa^ 

One can verify by a direct substitution that equation (JHJ) has the first integral 

ry Q Old/ A. . . 

11} = n A -n -n -, (10) 

? 1 + a 2 1 + a 2 1 + a 2 ' V ; 

where C is an arbitrary constant. Equation p()|l has the well-known solution in terms of 
elliptic functions. To write it down, we denote the zeroes of the polynomial in the right hand 
side of Eq. ffTUj) as ViiViiVz-, so that 

n\ = (n-pi){p 2 -n)(p 3 -n), Pi<p 2 <P3, (H) 

and suppose that n = p\ at £ = 0. As a result, we obtain 

n = pi + (p 2 - pi)sn 2 - Pi £/2; m) , (12) 

where sn(0; m) is the Jacobi elliptic sine and 

m = (13) 

Ps - Pi 

is the modulus. The constants A, B, C are connected with the zeroes Pi,P2,P3 by the relations 



2B 2C A 2 

Pl+P2+P3 = 7- J' PlP2 +P1P3 +P2P3 = -T- J, P1P2P3 = T- J" 

1 + a 2 1 + a 2 1 + a z 
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It is worth noticing that the components (u, v) of the velocity field are not determined 
unambiguously by the constants Pi,P2,P3 and a. Indeed, if n is known, we have only one 
equation (J7J) for calculation of u and v. Another equation can be added, if we restrict 
ourselves to the consideration of potential flows by imposing the condition 

Uy = V X , (15) 

which is consistent with the system (0) (see [T7]). and for a single-phase wave (0) yields at 
once 

au + v = D. (16) 

Here D is an additional integral of 'motion' which owes its existence to the Bernoulli theorem 
for the system (j5J). Equation (|16|) implies that the same spatial periodic profile of the density 
()12J) with the slope a can be supported by different potential velocity fields. If we fix u and v 
at some point, then the constant D becomes determined, as well as the velocity components 
everywhere. Just this situation occurs in the case of the soliton solution J7| where the 
velocity components are supposed to be known at |x| — > oo. Indeed, considering m = 1 in 
(|12j) and assuming n = 1, u = M = constant, where M is the Mach number and v = as 



\x\ 



oo we arrive at the oblique dark soliton solution obtained in |17j . 



q 



n = l z ; (17) 

cosh 2 [a/1 — q (x — ay) /(l + a 2 )] ' 

where q = M 2 /(l + a 2 ) and the velocity components are given by 

M(l + a 2 n) aM(l-n) , , 

(l + a 2 )n (l + a 2 )n v 1 

We note that we are not concerned here with the stability of the obtained periodic solutions 
which calls for further nonlinear modulation analysis. The stability of the oblique dark 
solitons (|17|) for M > 1 was established in ^7] numerically. 

It is important now to investigate the behaviour of the obtained periodic solution for the 
parameter values corresponding to some physically interesting asymptotic reductions of the 
system (JEJ). 



3 Small-amplitude nonlinear periodic waves 

As was indicated in J7| , if we consider a flow of BEC corresponding to small deviations from 
a uniform and homogeneous supersonic flow with n = 1, u = M, v = (M > 1), and make 
asymptotic expansions 

n = 1 + en\ + e 2 n 2 + . . . , u = M + eu 2 + e 2 n 2 + . . . , v — ev\ + e 2 v 2 + ■ ■ ■ , (19) 

where e 1 is a small parameter, then substitution of (fT9"j) into (0) followed by introduction 
of the scaled variables 

C = e 1 / 2 (x-Vy), r = e 3 / 2 y, (20) 
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leads, according to the standard reductive perturbation method, to relations 

= Vl = Ji ni ' V = VW^T, (21) 

where n\ obeys the KdV equation 

3M 2 M 4 . , 

riir ninic-\ rii ccc = 0- 22 

' 2VM 2 - 1 ,C 8VM 2 - 1 ,CCC V ; 

Its periodic solution is well known (see, e.g. [E]) so that the density profile can be expressed 
after returning to original (x, y)-coordinates as 



n 



2 

where 



1 - \m 2 e{\z - Ai - A 2 ) + M 2 e(A 3 - A 2 )sn 2 [v % - X 1 e l/2 {x - ay),m] , (23) 



q = y^T- £(Al + A2 + A3)M2 



4v^l ' m= A^' (24) 
and Ai < A 2 < A3 are the parameters arising in the finite gap integration method of the 
KdV equation (see [IS]). 

It is instructive to establish a direct asymptotic correspondence between fully nonlinear 
periodic solution (|T2*]l of the GP equation and its small-amplitude KdV counterpart (|2*Hj) . 
For that, we first represent the arbitrary parameters Pi,P2,P3 i n the form of asymptotic 
expansions in the small parameter e 

Pi = l + ep ( p + ... , 2 = 1,2,3, (25) 

and then substitute the expansions (fT9|l . (|2T|l . (J25|) into the periodic solution |T2|) to obtain 
the same solution (J23|) but parameterised by pf~' so that comparison with ([23]) yields 

Pi = 1 - ^M 2 £ (A 3 - Ai - A 2 ), 

p 2 = l-^M 2 e(A 2 -A 1 -A 3 ), (26) 
Ps = 1 - ^M 2 £ (Ax - A 2 - A 3 ), 
Then, using analogous asymptotic expansion for a, 

a = a (0) + ea {1) + . . . , (27) 

the relation (JZj) and the last relation in 1)14)1 . we recover the asymptotic expression ([21)) for 
the slope. Inverse expressions for A^s in terms of pjS, following from (|2l))l are, 

Al " NPe ' A2 " M% ' As " M% • (28) 

The soliton solution of the KdV equation corresponds to m — 1 in ()23)1 . In terms of the 
original parameters pj, we have p 2 = P3 = 1 (see ([T^jl . (J2E])), which implies by ()2~%]) that 
A x = A 2 = A and A 3 = 0. Then 

Pi = 1 + M 2 eA (29) 
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and, since p\ < p 2 < Pz ( see (HU)), we have p\ < 1 and thus A < 0. Now, from (j221) we have 
for m = 1 the small-amplitude dark soliton profile 



-M 2 e\ 
cosh 2 [V— s\(x — ay)] 

and from its slope is 



n = 1 ^ p=; -, (30) 



eAM 



2 



a = VM 2 - 1 (31) 

2y/M 2 - 1 V ; 



Since A < 0, one can see from (j31)l that a > gim, where om = V M 2 — 1 , i.e. the shallow 
(KdV) dark solitons always lie within (and close to) the Mach cone. 
If we introduce the inverse half-width k of the soliton according to 



k = 2 v /z eA, (32) 

then ()3())1 assumes a more conventional form 

n _ 1 M 2 k 2 

4cosh 2 [/t(x — ay) i '2] ' 

which also follows directly from the oblique dark soliton solution (j!7|) in the small-amplitude 
limit [IT]. 

Thus, we have established an asymptotic correspondence between the stationary periodic 
two-dimensional solution of the GP equation characterised by four independent parameters 
and the familiar three-parameter cnoidal wave solution of the KdV equation. 



4 Periodic stationary waves with large slopes in highly 
supersonic flow 

If the flow is highly supersonic (M ^> 1), then the system (jSJ) can be reduced to the NLS 
equation [TOj- Indeed, after introduction of the new variables 

u = M + Ul + 0(1/M), T = x/2M, Y = y, (34) 

where u\ — > as \Y\ — * oo, we arrive, to leading order in M _1 , at the system 

\n T + (nv) Y = 0, 
i , , ( n Y u yy\ n (35) 

\urr + VUiy = 0. (36) 

The leading term of the highly supersonic expansion of the potentiality condition (|T5|) to- 
gether with Eq. (f3T)j) implies u\ = 0. The decoupled from (f3T)|) equations (|33|) represent the 
hydrodynamic form of the ID NLS equation 

i^/ T + ^YY -2|*| 2 ^ = (37) 
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for a complex field variable 

-Y 



* = V^exp (i J v(Y',t)dY'\ 



(3f 



Note that this NLS equation, unlike the exact (1 + 1) reduction (J2J) of the GP equation, rep- 
resents a (0+2) asymptotic approximation and, in addition, contains a completely different, 
compared to (J2J), set of independent variables. Periodic solution of the NLS equation (J37)l is 
well known (see, e.g., [15J and for the density n = |\I/| 2 can be written in the form 

n = v 1 + (v 2 - ^i)sn 2 y/u 3 - v x (y - ^h^j , m , (39) 

where 

\ i \ i \ \ (A2 — Ai)(A 4 — A 3 ) 
si = Ai + A 2 + A 3 + A 4 , m = — — — — (40) 

(A3 — Ai){A4 — A2) 

and V\ < u 2 < z/ 3 are expressed in terms of the parameters A x < A 2 < A3 < A4 as follows, 

v x = |(Ai- A 2 -A 3 + A 4 ) 2 , 

^2 = |(Ai - A 2 + A 3 - A 4 ) 2 , (41) 



^3 — K^ 1 + A 2 — A 3 — A4 



,2 



and m = (u 2 — v\)j (z/ 3 - v\). 

For M ^> 1 the stationary wave has a large slope with respect to the y-axis, that is a 3> 1 
and the general solution (|12j) assumes the form 

n = pi + (p 2 -pi)sn 2 (y/p 3 -pi(y - x/a);m) . (42) 



Since the asymptotic solution (J42j) of the GP equation is characterised by four parameters 
Pi,P2,P3, a, the correspondence with the four-parameter (yi, v 2 , z/ 3 , si) periodic solution (JHnjl 
is readily established by a direct comparison: 

2M . . 

Pi = ^1, P2 = v 2 , p 3 = us, a= . (43) 

Sl 

The dark soliton reduction of the NLS cnoidal wave solution is obtained by putting 
A2 = A 3 , i.e. m = 1, in (J3~9"j) . To get a unit pedestal for the soliton it is convenient to choose 

Ai = -1, A 4 = 1. (44) 

Then, denoting A 2 = A 3 = A we get from (j4T7|) 

M , . 

a = — . (45) 

Thus, in the soliton limit we obtain 

1 — A 2 

n=l 2 —== (46) 

cosh [v 1 — \ 2 {y — x/a)] 
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which agrees with the asymptotic representation of a highly supersonic oblique GP soliton 
in ^7]. Without loss of generality we consider the waves in the upper half-plane (a > 0), 
then it follows from (|44j) that < A < 1. Thus, a > au and, therefore, the dark solitons in 
highly supersonic flows are always generated within the Mach cone. 

If the NLS soliton is shallow, then it must be consistent with the corresponding asymp- 
totic as M > 1 of the KdV soliton. Indeed, we introduce A 2 = 1 — M 2 k 2 /4, where Mk <C 1, 
and obtain the KdV soliton solution (|3*3*|) but now with a equal to 

a = M/y/l - M 2 k 2 /A = M + |M 3 k 2 (47) 

which is the approximation of Eq. (J24j) with chosen values of A, and M ^> 1. 

The connection between the soliton parameters A for the KdV and the "shallow" NLS 
limits is given by the expression 

*nls = 1 + M 2 s\ KdV . (48) 
Since Xxdv < 0, we have Xnls < 1 but 1 — Xnls ^ 1- 



5 Linear waves 



At last, let us consider the limit of linear waves propagating on a constant background (this 
case is more general than the linear wave limit within the KdV approximation as it does 
not imply the long-wave scaling (|2U|)). In this case Eq. (fP2|) with m <C 1 (i.e. p 2 — p\ <C 1) 
can be transformed to 



n 



1 - HP2 ~Pi) cos 



(49) 



where we have assumed that the mean background density is equal to unity, (pi +p2)/2 = 1. 
Further, in the limit of vanishing amplitude pi p\ and large wavelength we take 

Pi = p 2 = 1, p 3 = 1 + M 2 e, 



and obtain for M 2 e 1 



Vm 2 - 1 - 



M 2 e 



(50) 



(51) 



2VM 2 - 1 

Now we have a < cim, which means that the linear long waves, in contrast to solitons, are 
always generated outside the Mach cone. If we denote 



then we obtain the stationary linear wave solution in the form 

M 2 K 2 



Q, 



n 



Q cos 



x 



Vm 2 



1 - 



(52) 



(53) 



8\/M 2 - 1 , 

which maps to the tail of the oblique soliton solution in the KdV limit (|3*3|) by the change 
k i — > iK. 

As was noted above, the general formula (|49J) describes linear waves of an arbitrary 
wavelength and it can be mapped to the tails of the general soliton solution (fTTj) . 
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6 Conclusions 



We have obtained the family of exact fully nonlinear stationary periodic solutions of the 
2D Gross-Pitaevskii equation and studied in detail their particular asymptotic reductions 
corresponding to the solutions of the KdV and NLS equations. 

The obtained solutions provide a basis for further studies connected with the description 
of dispersive shock waves observed in recent experiments j2] of the flow of a BEC past 
obstacles as well as in numerical simulations ^7j. Some straightforward implications about 
the characteristic features of the wave patterns arising in the flow of a BEC past an obstacle 
have been made from the expressions for the slope a in the obtained asymptotic reductions of 
the full periodic solution. The deep solitons asymptotically described by the NLS equation 
have large slopes, while the shallow solitons obey the KdV equation and have slopes close 
the Mach cone, clm — 

VM 2 - 1, and they are located inside the Mach cone. The linear wave 
packets are always located outside the Mach cone. Detailed theories of these waves patterns 
generated by the flow of a BEC past obstacles will be developed elsewhere. 
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